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Abstract 

A strong-coupling expansion for the antiferromagnetic phase of the Hub- 
bard model is derived in the framework of the slave-boson mean-field approx- 
imation. The expansion can be obtained in terms of moments of the density 
of states of freely hopping electrons on a lattice, which in turn are obtained 
for hypercubic lattices in arbitrary dimension. The expansion is given for the 
case of half-filling and for the energy up to fifth order in the ratio of hopping 
integral t over on-site interaction U, but can straightforwardly be general- 
ized to the non-half- filled case and be extended to higher orders in t/XJ. For 
the energy the expansion is found to have an accuracy of better than 1% for 
U/t > 8. A comparison is given with an earlier perturbation expansion based 
on the Linear Spin Wave approximation and with a similar expansion based 
on the Hartree-Fock approximation. The case of an infinite number of spatial 
dimensions is discussed. 
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1 Introduction 



If the physics of strongly correlated fermions is to be described by the Hubbard 
model [ED], electronic structure calculations for the cuprates indicate that the rele- 
vant parameter regime is where the on-site interaction U is comparable to or larger 
than the bandwidth W of freely hopping electrons [[J . A natural strategy to try to 
describe this regime is to approach it from the two limiting cases of weak-coupling 
(U <C W) and strong-coupling (U ^> W). The simplest mean-field approximation 
for the Hubbard model, the Hartree-Fock approximation (HFA), is at first sight a 
weak-coupling approximation although, with some care, also a reasonable descrip- 
tion of the strong-coupling regime can be given Strong-coupling approaches 
have mainly been devoted to the one- dimensional case (see e.g. 0) and often take 
the limit U — > oo ||. A somewhat more sophisticated mean- field approximation, 
namely that based on the slave-boson formulation due to Kotliar and Ruckenstein 
0, is in principle not restricted to weak or strong coupling. In fact, it was shown 
that this so-called slave-boson mean-field (SBMF) approximation is equivalent to the 
Gutzwiller approximation to the Gutzwiller wave function fpTOf] , the latter of which 
improves upon the HFA especially in the intermediate-coupling regime. This SBMF 
approximation then allows to approach the interesting intermediate-coupling regime 
from the strong-coupling side. Such an approximate strong-coupling approach can 
be particularly helpful since more rigorous Quantum Monte Carlo calculations be- 
come increasingly cumbersome for stronger coupling [II]. 



In a previous work [12], we compared the HFA and SBMF approximation for 



various simple magnetic phases, as well as computed the effective hopping and spin 
stiffness for the two-dimensional (one-band) Hubbard model. In this paper, I derive 
a large-?/ expansion within the SBMF approximation for the Hubbard model in 
arbitrary dimension. Here, I restrict myself to the antiferromagnetic phase at half- 
filling, but the extension to doping the antiferromagnet with electrons or holes is 
straightforward (though tedious). 

In Section 2, the slave-boson mean-field approximation is briefly introduced and 
the resulting set of integral equations for the antiferromagnetic phase is given. A 
large- expansion is derived in Section 3 in terms of moments of the density of 
states for freely hopping electrons on a lattice, and compared to earlier numerical 
and analytical work in the literature. In Section 4, the expansion is made explicit 
for hypercubic lattices of arbitrary dimension. Finally, I discuss the convergence of 
the series and possible extensions. The Appendix contains useful formulae for the 
moments of the above mentioned density of states. 



2 Slave-Boson Mean-Field Approximation 

The Hubbard Hamiltonian is pj: 

H h = -tY, c]+s,*Cj* + U J2 n rt n ih (!) 

3,8," i 

where Cj a is the annihilation operator for an electron at site j with spin a. Neigh- 
boring sites of site j are denoted by j + 5. t is the hopping integral, U the on-site 
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interaction between electrons of opposite spin, and rij a = c^ a Cj a is the occupation 
number operator. The slave-boson approach of Kotliar and Ruckenstein consists of 
introducing bosons for each of the (four) possible electron occupancies of a site. The 
electron creation and annihilation operators are then modified such that always the 
(one) boson corresponding to the electron occupancy is present at each site. In a 
functional-integral description of the new Hamiltonian the slave-boson mean-field 
(SBMF) approximation is when the Bose fields are independent of (imaginary-) 
time. Different phases can be considered by assuming different forms of position 
dependence of the Bose fields. For more details see Refs.|| |12|]. The SBMF approx- 



imation is an improvement over the Hartree-Fock approximation since it takes some 
local correlations into account; in particular, the density of doubly occupied sites is 
an independent parameter to be optimized. 

The set of equations describing the antiferromagnetic (AF) phase within SBMF 
is given by [p~2]j : 



m, = 2A / As M(C 



J » (e 2 + A 2 ) 
q s H ( £ 2 + A 2^ 



[q*]d f 00 ^ M(e)e 
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poo 

n = 2 deX(e) , (5) 
J n 

where A is an internal ( "renormalized gap") parameter; p, an effective chemical 
potential; and the band renormalization q s is a function of density n, sublattice 
magnetization m s , and d (d 2 is the density of doubly occupied sites), and given by: 

q s (n, m B , d) = z(n, m s , d)z(n, -m s , d) , (6) 

where 



[1 -n + d 2 )(n + m s -2d 2 ) +dV n — in* — 2d 2 
z(n, m s , d) = 1 . (7) 

sj{n + m s ) (l - s^) 

[q s ] a is the derivative of g s with respect to a (explicit expressions can be found 
in Ref . [fL2]| ) and Af(e) is the density of states (DOS) of freely hopping electrons. 
The dimension- dependence resides solely in M{e) and since this DOS applies to a 
non-interacting system, I will be able to derive results for arbitrary dimension in 
Section 4 (see also the appendix). Note that the above set of equations reduces 
to the familiar antiferromagnetic or spin-density wave (SDW) solution in the HFA 
if g s = 1 || . In that case, one has one self-consistency equation (0), which is the 
gap equation with A as gap parameter A (= Um s ), instead of three self-consistency 
equations ©-©. 
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The energy of the AF state (per site) and the spin stiffness p s are given by: 



e AF = -2q s [ &eN{e)\[e 2 + A 2 + Ud 2 + Xq s m s 
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00 . A/YeV 
as 



e 2 + A 2 



1/2 



zlz 2 



A' 



e 2 + A 2 



3/2 



(8) 
(9) 



Af v (e) is the weighted DOS, which, as well as the combination z+z~, is specified 
in Ref . ||1 2|| . I will not repeat the expressions here since I will now restrict to half- 
filling (/2 = 0, n — 1) for which case z + z_ equals zero. For half-filling, the band- 
renormalization reduces to: 



q s (l,m s ,d) 



2d 2 



VI - 2d 2 + m s + VI - 2d 2 - m s 



(10) 



3 Large-?7 expansion in terms of moments 

The above set of equations (0H([5|) can be solved numerically if the DOS for freely 
hopping electrons is given. For instance, in two dimensions N(e) is known analyti- 
cally and numerical solutions to (0)— were obtained in Ref. [|T2H ■ Now I will derive 
solutions to (|2|)-(|5|) for n = 1 in the form of a systematic series expansion in 1/U. 
This expansion is derived in two steps: (i) From the explicit form of the equations 
above (or from numerical solutions) it can be seen that large U corresponds to large 
A. Therefore, I first make a large-A expansion for all quantities of interest. Sub- 
sequently, (ii) the A-expansion for U (following from equation (|])) is inverted and 
substituted in the large-A expansion for all the relevant quantities obtained before. 
In this way, one obtains the desired large-?/ expansion. 

By expanding the square root in the denominator of the integrals in (0)-© for 
A ^> £ one obtains a large-A expansion in terms of moments M n of the DOS: 



M„ 



dsAf( 



EE 



'IV. 



After some work one arrives at the following expansions in q = 1/A 



d 



1 - \M 2 q 2 + |M 4 g 4 - ^M 6 q 6 + ^M 8 g 8 + 0(q 



V0\ 



4Mf + 3M 4 ., 



q° + 



fW 2 
16 



M 2 + f M 4 - 



\^jM 2 q~ 

-1 

2048M 2 5/2 

736M 2 3 M 6 - 120M 2 M 4 M 6 + 280M 2 2 M 8 j q 7 + 0(q 



9Mj 
16M 9 2 



+ 



192M 2 6 + 368M 2 4 M 4 + 364M?M 4 2 + 27M 4 3 + 



M, 



2 „2 



M| - 7M 4 
16 



9M| - 10M 2 M 4 + 12M 2 /M 2 + 26M 6 



64 



g 6 + 0(q 8 



5M 6 
2M, 



(12) 



(13) 



(14) 
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u 



eAF 



+ f M 2 q + 



Ps 



Q 

-\M 2 q + 
M 2 



11M 2 2 



21M 4 



3M| + M 4 



g 3 + 



M| + 2.U; 

32 



g 3 + 



g 3 + 0(g 5 ) , 

"M| - 22M 2 M 4 - 2M 6 
32 



-M| + 9M 2 M 4 + 6M 6 



(15) 

q 5 + 0(q 7 ) ,(16) 



128 



9M 2 4 - 12Mf M 4 + 26M| + 32M 2 M 6 + 20M 8 
512 



g 7 + 0(q 9 



(17) 



For notational convenience the hopping integral t is taken equal to 1 in this section. 
For m s , d, and q s I have listed precisely the number of terms that is needed in order 
to obtain the number of terms (three) that is listed for U. For caf and p s the terms 
are given that can be found using the listed terms for m s , d, q B , and U. Note that 
this implies that, to obtain the leading order in q for U, one already needs three 
terms for m s and two terms for d. 

Inverting the expansion for U gives: 

q = 2p + 10M 2 p 3 + (122M 2 2 - 42M 4 )p 5 + 0(p 7 ) , (18) 

where I have introduced the notation: p = 1/U. Substituting this into the q- 
expansions, one obtains the following large-[7 expansions: 



d 2 
Qs 

eAF 

Ps 



1 - 2M 2 p 2 



20M| - 6M 4 



p 



294M 3 - 204M 2 M 4 + 20M 6 J p 6 + 0(p 8 ) , 
M 2 p 2 + \QMl - 3M 4 1 p 4 + [75M 3 - 70.U..U; + 10M 6 ] p e + 0(p 8 



1-M 2 p 2 - 



11M 2 - 7 M 4 



P 



176M 3 - 192M 2 M 4 + 12M 4 2 /M 2 + 26M 6 p 6 + 0(p 8 



-M 2 p 
M 2 



2.1/.7 - M d 



■2 M; - M; 



p 



15M 2 3 - 14M 2 M 4 + 2M 6 p 5 + 0{p 7 



(19) 
,(20) 

(21) 
(22) 



p 3 + § 15M 2 3 - 14M 2 M 4 + 2M 6 p 5 + C(p 7 ).(23) 



A consistency check on the above results is that the following general identity still 
holds for the expansions above: 



de 



AF 



dU 



d 2 



(24) 



Because in the SBMF at half-filling the spin stiffness is given by — | times the average 
kinetic energy [l^j, we also have: 



Ps 



I (e AF - Ud< 



(25) 



Combining the latter formula with (|4]) it immediately follows that the expansion 
for p s follows directly from the expansion for eAF, as is also seen in the way the 
expansions are written above. 
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It is of some interest to compare the present result with results from the HFA 
and other expressions in the literature. It is far easier to derive the corresponding 
expression to (22) in the HFA since there is only one consistency equation (the gap 
equation) instead of three, and there is no band-renormalization function q s which 
requires perturbative expanding. Following the same procedures as for SBMF, I find 
for the HFA: 



e AF = -M 2 p - M 2 2 - M 4 



V 



AM 3 



6M 2 M 4 + 2M 6 p 5 + (D(p 



(26) 



Comparing with (p2|), one sees that HFA and SBMF give the same coefficient for 
the leading order contribution, but that there are differences for the higher order 
contributions. Note that also for the higher order contributions the term in the 
coefficient involving the highest moment agrees between HFA and SBMF. Since the 
moments are always positive, it is also clear that for large U the energy in the SBMF 
approximation is always lower than in the HFA. 

In dimensions D = 2 and D = 3 (square and simple cubic lattices), the results can 



also be compared with a large-?/ expansion by Takahashi . He derives a rigorous 
expansion (i.e. without making approximations) for the half-filled Hubbard model in 
terms of spin correlation functions of the s = § Heisenberg model. These correlation 
functions are then evaluated using the Linear Spin Wave (LSW) approximation. 
Below the results of the latter approach are compared to those of HFA and SBMF 
for D = 2 and D = 3 (the necessary moments are easily evaluated and may for 
instance be found in the appendix): 



HFA 
SBMF 
LSW 



D 

Ap + 20p 3 



2 

192p 5 + 



-6p 



D = 3 

54p 3 - 1344p 5 + 



- Ap + Ap 3 + 256p 5 + ■ 
4.63p + 34.6p 3 + • • • 



-6p + 18p 3 + 600p 5 + 
-6.58j9 + 65.6p 3 + --- 



(27) 



It is clear that for large enough U the LSW expansion will give the lower energy. 
The leading order coefficient is proportional to the ground-state energy of the s = 1 
Heisenberg antiferromagnet, which is known to be very well approximated by LSW 
15L [16|]. Both the HFA and SBMF do not do better than the mean- field result for 



the leading order coefficient. In view of the fact that SBMF is an improvement over 
the HFA, it is somewhat remarkable that the HFA results for the next-to-leading 
order resemble the LSW results more than the SBMF results do. Also the different 
sign of the coefficient of p 5 between HFA and SBMF is notable. Further discussion 
of the range of validity of the expansion is given below. 

To end this section a quantitative discussion of the large-?/ expansion is also 
given. In Table I, the ground-state energy of the AF phase at half-filling as ob- 
tained from the large-?/ expansion within SBMF (e^p xp ) is compared to the "exact 
SBMF" result [e^p MF , found by solving (^)-([5]) numerically] and to variational Monte 
Carlo results using an (antiferromagnetic) Gutzwiller wave function (e^p YMC , from 
Ref.0) for D = 2. For D = 3, only e^ p and e^ VMC are compared. Note that 
the large-?/ expansion approximates the "exact SBMF" result very well for U/t > 7; 
it is expected that similar agreement is found for D = 3, where it is more involved 



6 



to compute e AF MF since Af{s) is not a known analytical function as it is for D 



TS | . This is consistent with the fact that from the coefficients in the explicit expan- 
sions for D = 2 and D = 3, formula (|27j), one would estimate the series to converge 
well for U/t > 8 and U/t > 6 for two and three dimensions, respectively. SBMF 
is equivalent to the Gutzwiller approximation to the Gutzwiller wave function || 
and this approximation becomes exact in the limit of an infinite number of spatial 
dimensions [19. 20 1. From the comparison with the GWVMC results it can be seen 
that the approximation is quite good already for D = 2. In view of the facts, that 
for D = 3 a very small lattice was used in the GWVMC calculations and a larger 
lattice will even raise the energy somewhat (the effect will be larger for smaller U), 
the agreement for D = 3 can even be called excellent. So the clear advantage of 
the large-?/ expansion is that it can give results that are, for U/t>8, within 1% of 
the "exact SBMF" and GWVMC calculations (which especially for D = 3 are much 
more involved and computationally demanding) by means of the simple formula (f2^). 







D = 2 




D 


= 3 


u 


„SBMF 
— e AF 


Uexp 
— e AF 


„ GWVMC 
— e AF 


Uexp 
— e AF 


„ GWVMC 
— e AF 


6 


0.623639 


0.615226 


0.629(3) 


0.8395 


0.886(5) 


8 


0.485104 


0.484375 


0.493(3) 


0.6965 


0.704(7) 


10 


0.393528 


0.393440 


0.401(4) 


0.5760 


0.579(6) 


12 


0.330002 


0.329990 


0.336(5) 


0.4872 


0.491(7) 


16 


0.248782 


0.248779 




0.3700 




20 


0.199420 


0.199420 




0.2976 





Table I. Numerical results in D = 2 and D = 3 for the ground-state energy of the 
antiferromagnetic phase at half-filling. e^p xp is evaluated using the large-f/ expansion 
flS^) and e|f T MF is obtained from <g) after solving <g)-(g). The GWVMC results are 
taken from Ref.jlj]] and were obtained on 20 x 20 and 6x6x6 lattices for D = 2 
and D — 3, respectively. Between brackets the statistical error in the last digit is 
given. 



4 Large- U expansion for arbitrary dimension 

In the appendix, it is shown that the moments M n can be obtained as a function of 
dimension D [formula ( |A.9| )]. If one substitutes these expressions in Q2"2| ) and ( [26] ) 
one obtains a simultaneous l/U and 1/D expansion for the energy of the AF phase: 

2Dt 2 ( 3 \ ADH 4 ( 24 20 \ 8DH 6 

SBMF e AF 



HFA e AF 



2Dt 2 


+ 




ADH 4 


U 


U 3 


2Dt 2 


+ 




ADH 4 


U 


U 3 



+ 



+ 



24 20 
3+ D 



-16 



36 
15 



20 
D 2 



U 5 
8D 3 t 6 



+ ••(23) 

+ •■(29) 
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where I have included the hopping parameter t again. The above formulas are well 
suited to discuss the limit of an infinite number of dimensions. This is of interest 
since, as noted above, SBMF is equivalent to the Gutzwiller approximation to the 
Gutzwiller wave function; an approximation which becomes exact in the limit of 
an infinite number of spatial dimensions. Thus the present expansion allows us to 
scrutinize the approximation explicitly for lower dimension. In order for the limit 
D — > oo to be meaningful, one has to introduce a scaled hopping parameter t* = 



ty2D [19||. In terms of this scaled parameter the above SBMF large-f/ expansion 



assumes the form: 



!£ = _!lJ,_lvn 3 J- 3 4- «) ( £) 5 + . . . ,30) 



t* U V 2D J \UJ V D D 2 J \U 
In Table II, I evaluate (SO) for D = 2, 3, 4, and oo. 



u* 


-4fP = 2) 


-e* AF (D = 3) 


-elpp = 4) 


-e* AF (D = oo) 


2 


0.34375 


0.35069 


0.36719 


0.46875 


4 


0.24219 


0.23947 


0.23853 


0.23730 


8 


0.12439 


0.12394 


0.12373 


0.12314 


12 


0.08317 


0.08303 


0.08296 


0.08277 


16 


0.06244 


0.06238 


0.06235 


0.06226 


20 


0.04997 


0.04994 


0.04992 


0.04988 



Table II. Scaled ground-state energy of the AF phase at half-filling e^ F = e^p/t* 
within SBMF as a function of scaled interaction U* = U/t* for various dimensions 
D (with t* = tV2D). 



From Table II, one sees that for scaled interaction U* > 4 (which corresponds 
to U/t = 8 for D = 2) already the result for D = 2 and the infinite-dimensional 
result coincide remarkably well. The simultaneous large-f/, large- .D expansion shows 
in detail how well the Gutzwiller approximation, and therefore SBMF, reproduces 
the Gutzwiller wave function for the various dimensions. It would have been in- 
teresting to compare the coefficients in the present results with those obtained by 
Gebhard in his 1/D-expansion for the half- filled Hubbard model [ 20[ . Unfortunately, 
a direct comparison is not meaningful, since Gebhard only considers the paramag- 
netic Gutzwiller wave function and does not treat the antiferromagnetic phase at 
half-filling in a 1/D-expansion. 
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5 Discussion and Conclusions 



There is no fundamental reason why the large-?/ expansion given here cannot be 
extended to higher orders in t/U . With a little more effort (especially the coefficients 
in the large- A expansion of d in formula (|H|) require some care) more terms can be 
obtained. However, as already remarked in Section 3, the present series reproduces 
well the "exact SBMF" result for U/t > 7 in two dimensions, and this agreement is 
expected to be even better in three dimensions. Also, for smaller values of U/t the 
contribution from the {t/U) 5 term becomes larger than that from the (t/U) 3 term, 
which is an indication that the intrinsic radius of convergence of the series is close 
to U/t = 7. Therefore, extending the series can not be considered very useful. 

A more interesting extension of the present series is to go off half-filling, i.e. 
n^l. Again there is no fundamental reason preventing this, but now the task is 
quite a bit more laborious for a number of reasons. Not only are we dealing with 
the case p, ^ and consequently partial moments are required [see (II)], which can 
only be computed numerically, also the expression for the band renormalization q s , 
formulae (|6])-([7|), needs to be expanded for n ^ 1, which implies that non-rational 
coefficients in the expansion in terms of (partial) moments appear. Furthermore, in 
formula ([5|) for p s also the term with the combination z + Z- needs to be taken into 
account and perturbatively expanded. If similar large-?/ expansions are possible for 
more complicated phases, like spirals or domain walls, this would allow for a study 
of aspects of the phase diagram off half-filling. 

In conclusion, a large-f/ expansion is derived for the antiferromagnetic phase of 
the Hubbard model within the framework of the slave-boson mean-field approxima- 
tion. Even though such an expansion constitutes in a sense an approximation of an 
approximation, the resulting analytic expression is capable of reproducing very well 
results of elaborate Monte Carlo calculations. Furthermore, the expansion allows 
for an explicit study of the limit of a large number of spatial dimensions, since the 
occurring moments of the non-interacting density of states are obtained for arbitrary 
dimension. 
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Appendix 



In this appendix some general expressions are derived for the moments M n of the 
density of states M{e) of freely hopping electrons on hypercubic lattices: 



deJ\f(e)e n = ^ ; / ■ • • / dkx . . . d£; D [cos fci H h cos k D ] n . 

-co {27l) 1J J-tt J-tt 

(A.l) 

Although such moments are frequently used and calculated in the literature I have 
never encountered the general analytic expressions and simple formulae to be given 
below. First, I will derive simple expressions for all moments in the lower dimensions 
D = 1,2, and 3. Then I will derive expressions for the lower moments in arbitrary 
dimension D. 

By employing the multinomial generalization of the binomial it is possible to 
further work out ([A. 1|) : 

n , n) 

M n = (2t) n J2 —— :q m ...q nD , (A.2) 

n lr ~^=o n x \...n D \ 

where the prime on the sum denotes the restriction: n\ + • • • + ud = n and 

q m = - dkcos m (k) = \ m ^ (A.3) 

2n Jo { m odd 

From this general form it is clear that all moments with n odd are zero for any 
dimension. All formulae below will be for moments with n even. For D — 1 the 
restriction allows for only one term in the sum and the result is trivial: 

(n - IV 

M n = (2t) n [ D = 1 . (A.4) 

A more remarkable result is that the moments for D = 2 are exactly the squares of 
the moments for D — 1 (if t is put equal to 1): 



M n = f 



(n-DW 1 ' 2 



nil 



D = 2 . (A.5) 



The proof of (|A.5|) follows directly from ( |A.2|) by some manipulations with the 
factorials and employing the identity 



N 'N\ 2 (2N 



sw-uj- (A - 6) 

For D = 3 I have not been able to obtain an explicit expression for M n , but the triple 
sum in ( |A.2j ) can be reduced to one simple, unrestricted sum. Using the relations: 
(m — l)!!/m! = l/m\\ = 2~ m ^ 2 / \m/2)\ (m even), one can rewrite (|A.2|) (for n even) 
as: 

n/2 1 

M ^ n \^Wlm^' (A ' 7) 
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where the prime indicates the constraint: u + v + w = n/2. This expression was 
already given by Brinkman and Rice for the moments of the density of states of 
one hole in a ferromagnetic background |[22|| . It can however be simply reduced to 
a single sum by incorporating the constraint and a subsequent application of (|A.6| ). 
The result is: 

n/2 (n-2uV 
M n = n\y 1 h j D = 3. (A.8) 



Formulas (|A.5|) and (|A.8| ) make for much easier evaluation of the moments than is 



usually found in the literature: e.g. in Ref.p3| moments up to n = 22 and n — 16 



are computed for D = 2 and D = 3, respectively, using the method of counting the 
number of paths that return to the origin after n steps. I note that also finding an 
explicit expression for D = 3 may be possible on account of the (empirical) fact that 
the moments for D = 3 always contain the moments for D = 1 as a factor. 

The general expression ( |A.2[ ) shows that finding the terms contributing to M n 
amounts to finding all the combinations of D even numbers (including zero) that 
add up to n. For the lower moments the number of possibilities is not large and 
can be expressed in terms of D using combinatorial arguments. I have found the 
following results: 

M 2 = 2D 

M 4 = QD{2D - 1) 

M 6 = 20D(6D 2 - 9D + 4) (A.9) 

M 8 = 7(LD(24D 3 - 72D 2 + 82D - 33) 

M 10 = 252D(14400D 4 - 143400D 3 + 501050D 2 - 715225D + 343176) 
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